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The confinement of electromagnetic field is studied in axial symmetrical, warped, 6D World Brane, 
using a recently proposed topological abelian string vortex solution as background. It was found, 
that the massless gauge field fluctuations follow 4D Maxwell equations in the Lorenz gauge. The 
massless zero mode is localized when the thickness of the string- vortex is less than 5/3/47re^«^ there 
are not others localized massless modes. There is also an infinite of non localized massive Fourier 
modes, that follow 4 dimensional Proca equations with a continuous spectrum. To compute the 
corrections to the Coulomb potential, a radial cutoff was introduced, in order to achieve a discrete 
mass spectrum. As main result, a correction was found for the 4D effective Coulomb law, 
the result is in correspondence with the observed behavior of the Coulomb potential at nowadays 
measurable distances. 
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'sj- . I. INTRODUCTION 
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In the brane world view, the universe is considered a 3-brane contained in some topological defect. This 3-brane is 
embedded in non compact, large, extra dimensional spacetime with finite volume 111 usually with warping metric [2j| in 
the directions of the additional dimensions. Most of the research had been done [3| on domain walls in 5 dimensions, 
Qh] due to the fact that this is the simplest kind of topological defect and exact solutions are available for both the metric 
^ • and scalar fields. Important advances had been made on domain wall's brane world, as possible explanation to the 
hierarchy problem Q , localization of gravity [1] , and the confinement of some fundamental fields Q as chiral fermions 
, but there are still several unsolved problems as the localization of gauge fields Q , stability under casimir energy 
' or quantum fluctuations [6|] and the lack of a complete supersymmetric version. 
^ I Its is widely known that in 5 dimensional brane worlds, the gauge and spinor fields can scape from the 3 brane 
universe into the bulk. Yukawa coupling gives a natural way of localize fermionic fields on topological defects [3| 
but the confinement of spin 1 fields remains as an open problem. Some mechanisms as DS 3 and DGS of quasi 
localization has been proposed. These mechanism seem to unnaturally force, or increase the effect of the "wall" over 
the volumetric "bulk" terms. To do that, they introduce two different coupling constant for the electromagnetic field 
and very strict restrictions on the fields and parameters. Some recent calculations indicate very strict bounds on the 
■ wall thickness of many orders of magnitude thinner than the Planck scale 10] to achieve a 4 dimensional behavior. 
y—{ , The loss of electromagnetic radiation at low frequency and the non preservation of charge are also drawbacks related 
1" ' to the confinement of gauge fields in domain wall brane worlds. Although this drawbacks can be controlled in several 
. ^ ' ways, they are also unwanted features for a theory modeling the universe. 

Some years ago, brane worlds had been proposed over more complicated topological defect than domain walls. For 
?H ■ example world brane had been proposed for the 6 dimensional vortex pT|[T2j and 7 dimensional monopole [Tsjfl^. 
The confinement of spin 1 gauge fields [isj seems to be possible to achieve in strin g-vo rtices in 6 dimensions: it has 
been found that electromagnetic gauge field fluctuations has a localized zero mode ]l^ [11] , however the graviphoton 
zero modes (spin 1 fluctuations of the metric) are not normalizable due to its behavior either at infinity or very near 
to the string- vortex [17] 18]. These works are based on numerical or asymptotical approximations, because there are 
not known general analytic solutions for curved 6D string-vortices and very few exact solutions had been reported 
[2^. Moreover, almost nothing is known about the non zero modes and its influence on the propagators, potentials 
and interactions. 

Recently, new Randall Sundrum scenarios in 6 dimensional curved space time based on exact topological solutions to 
the abelian higgs model [1^ had been reported. These solutions correspond to different Einstein Maxwell field vacua, 
with scalar kink solitons and auto-interaction potential, that are in fact, electromagnetic uncharged string-vortices 
with non trivial winding number. These topological abelian strings exhibits a localized spin 2 zero mode, that give s 
rise to the newtonian gravity, while the rest of the modes gives account of a correction to the potential [Tl| . 
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Is the purpose of these paper to study the confinement of spin 1 gauge field, on a background of a 6 dimensional 
topological abelian string vortex, in order to establish the localization or not, of all the photon's fourier modes and 
calculate the corrections to the Coulomb law. 

The organization of the paper is as follows: section II explains the model: the topological abelian string vortex. 
Section III reviews the confinement of the zero mode of linearized gravity and the non confinement of the other 
modes. In section IV, the gauge fixing conditions and the decoupling from the graviphotons and graviscalars field are 
discussed and the 4D Maxwell equations recovered in a perturbative approach using the topological abelian string 
as background. In section V, the gauge field equations are expanded in Fourier modes. A localized, normalizable, 

solution to the massless zero mode is obtained numerically in the thin string approximation: 6 < ^^^2^2 it is shown 
that there are not other localized massless modes. In section VI the massive modes are studied, a regularization is 
proposed based on the introduction of a radial cutoff and the corrections to the Coulomb law calculated. It is the 
main result of this paper a correction to the Coulomb potential. In section VII the range of validity of this 
corrections and its adjustment to known experimental bounds is discussed, finally some conclusions and remarks are 
presented. 



II. A TOPOLOGICAL EINSTEIN ABELIAN STRING- VORTEX SOLUTION. 



We start from the 6 dimensional action 



1 dx^V^ 




+ j dx^v^ 




. 2x 





-FabF^^ - V 



(1) 



where the first integral is the 6D Einstein Hilbert action, with a bulk cosmological constant A and Db — Ob ~ icAb 
is the covariant derivative. 

In what follows we will set y — t^txt — 1, where Mr is the 6D Planck mass and we will follow the notation in 

MM- 

Here we will look for a geometry 2+3+1 composed by a 3-brane that contains the 3+1 physical universe and 2 
extra dimensions, where we can choose coordinates (r, 6) and the metric is given by: 



M\r) i^^^dxi'dx" - Ri L^{r)de^ - dr\ 



(2) 



with ry^y = (+, — , — , — ). In this context r G [0, 00) and 9 G [0, 2tt] are the coordinates of the extra dimension, R^L{r) 
acts as a radial factor and Ad{r) is a warp factor, rapidly diminishing when moving away from a 4 dimensional 3-brane 
located at r = 0. 

Assuming that the scalar and gauge fields in ([T]) depends on the extra coordinates as in the Nielsen Olesen ansatz 



x^,r 



= V f{r) neZ (3) 
n — P(r\ 

A,(x^r,^) = (4) 

e 

Ar{x^\r,e) = 0. (5) 

A^{x^',r,e) = (6) 

where the w = 1 is dimensional factor with units lenght"^. Performing variations in the action ([T]) we get the curved 
version of Nielsen Olesen vortex equations [T2j : 

^ + (4m + 0^-^-^^ - 

^ + (4m-0^-^V/2p = 0, (8) 
ar^ ar 

where 

P = n [1 ~ a{r)]; m = ln[M(r)]; I = ^\n[L{r)]. 

String-vortex solutions are classically obtained with the boundary conditions: 
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/(r^0) = 0, f{r^^) = l, 

a{r — 0) = 0, a{r — > cxj) = 1, ^ ' 

by numerics and asymptotical methods. Very few solutions to the equations (j7l8p are known. Indeed only the 
Bogomoln'yi solution for flat space in the critical case [2^ is known to be exact. Recently in [l^[23| the 
boundary condition a{r) = 1 was used to find new exact solutions, so we will assume here: 

P(r) = 0, V re(0,cx)) (10) 

and the equation ([5]) is automatically satisfied. 

Boundary condition (I10[) lead to a apparent singularity at in the vectorial field because the direction ug is not 
defined at r=0, and the potential 1-form is locally defined by 

A = Agdx^ = Aed0. (11) 

The model ([1]) is invariant under the group U(l) of local gauge transformation: 

^ e^^(^)(j) Ab = Ab + -dsHx) 



e 



So, when the phase of the scalar field is choose as 



d{x,r,e) = n9 =^ SAg = - 

e 



that is exactly Q when a(r) = 1 and is valid for r 7^ 0. 

Equation ([Tl]) generates a null electromagnetic field Fab = 0, and although this type of string- vortex has neither 
electric nor magnetic charge, it is a topological vortex solution because it still has a non trivial integer winding number: 



I A, neZ (12) 



C 

where C is any closed curve around a "string" at r = 0. 

It is not possible to continuous pass or continuously deforms a curve with a particular value n to a curve with a 
different n. For each homology class of curves C we will get the same integer value for n. Each different value n will 
correspond to a class of curves C and labels an specific topological vacuum. 

This kind of string- vortex solutions will be referred as "topological string- vortices" [ISl and they are useful to 
construct 6 dimensional brane world. In order to obtain a feasible Randall Sundrum scenario in 6 D, we still have to 
prove that that the metric factor (r) is warped and the potential shows breaking of the symmetry. 

Einstein equations are obtained performing variations of the metric in the action and jointly with ([7]) conforms a 
system of 4 coupled non linear differential equations and 4 variables (m, /, /, V). With straightforward combinations, 
the complete system could be written as: 

f" + {4m + l)f-^-^V = 0, (13) 

I' + {4m + 1)1 = -^V, (14) 

m' + (4TO + /)m = (15) 

m' + -ml = -^{f'f- (16) 

Where the 6D (bulk) cosmological constant was absorbed into the redefinition of the potential V — !■ V + A and 
although X = 1 E^J^d v — 1, had been introduced in the equations to match the dimension of the system of units. 
Equations dH]) and p5)) implies m = I, so the system reduces to: 

r + Sm./'-l-^F = 0, (17) 



5m2 (18) 
m' = -1(/')^ (19) 
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for a given V{f) the system (|17I18I19|) is not longer independent, and always one of the equations can be obtained 
from the other two. But solving this system is not an easy task, due to its coupled non linear nature. 

We will follow here the approach developed in [l^ and [l^ , more details could be found in [1^ . Instead of try to 
solve for a given V{f), we will give a probe function f(r) that accomplishes the boundary conditions. 

So, from (IT9l) it is easy to obtain m as: 



kRS, 



(20) 



where k^s is an additional Randall Sundrum constant warp factor. Then, from equation (jl8l) the potential V{r) 
could be obtained as a function of r. In order to obtain the interaction potential V{f) we must solve r{f) from the 
probe function /(r) so: 



V{f) = V{r{f)), 



(21) 



Of course a solution to be physically acceptable must have a stable potential with spontaneous breaking of symmetry, 
as we will show immediately. 

Equations p7ll8ll9p are very similar to that of 5D domain walls in [IJl , so we will try: 



/ = /o arctan (sinh i3r/S) , /o = 2^5, 

using ([20]) we get 

m — — /3tanh { jSr / 5) , where we set kf>s = 0, 
That corresponds to the metric "warped" factor 

M{r) ^ L{r) ^coshr\l3r/5), J > 0, ^ > 0. 
From equation ([TSl) we obtain the potential V — V{r) as: 

Vir) = ^[(1 + h5)sech^{l3rl5) ~ 5(5], 





(22) 
(23) 
(24) 

(25) 



and solving r as a function of / by means of we have cos^( j-) — sech^{^). So finally the interaction potential 
for the scalar field is 



n/) = ^[(1 + 5(5) cos^ (///„) -5^]. 
ox 



(26) 



This is an stable potential, with spontaneous symmetry breaking minima, where the scalar field interpolates 
smoothly between / = ivrv^ in AdS space-time with cosmological constant A = — 10/3^/x as could be seen in 
FigUl 






FIG. 1: Graphic of the metric factor M, scalar field (left), and the energy density (right)vs r. Graphic of the auto-interaction 
potential vs / (center), /o = 1, 5 = 1/4, /3 = 0.1, X = 1 
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III. LOCALIZED GRAVITY ON THE TOPOLOGICAL ABELIAN STRING 

We will now study linearized spin 2 metric fluctuations from the metric (|24p given at first order by 

gyr ~ M'^[riyr + Kr] (27) 
we look for eigenfunctions with separated variables of the form; 

Ku {x^,r,e) = hr, (x^' ) 0^ (r)e''=^ , (28) 

K 

Then the gravitation equation from ([T]) splits into two equations as in [Tlj 

d^d^Krix") = ml Krix") (29) 

= (30) 



d 



M'^Ldr 



M^LiM^^icf,,) 



= - ^ (31) 



where is the constant of the separation of the equation for the gravity in variables (x") and {r,9). Here the 
eigenvalues /x are effective mass factors that depends on TOq and k. 
Looking for solutions of the form of flat waves: 

equation (|29p could be interpreted, as massive graviton with mass term PvP'^ = > in 4 dimension. Equation 
P0|) could be though, as the radial equation of a massive scalar mode, in curved 2 dimensional space (r, 6) with an 
effective mass given by /x^ > 0, that is also non negative because also paP^ > in 6 dimensions. 

Upon substitution of and using M{r) = L{r), we get the equation for massive metric fluctuations: 

- + 5/3 i!mh{Pr/5)(j)'^^ - cosh(^r/(5)<?!)^ = 0, (32) 

So integrating the equation (15^ when /x^ = 0, we obtain the massless mode 

<j)^^ko + ki [ drcos\i^\prl5). (33) 



with integrations constants ko and ki. 

As coa\^^ [j3r / 5] is monotonous growing, we must fix ki — 0, in order to render (po normalizable. That is consistent 
with the boundary conditions 

M^) = <t>'^^) = 0, (34) 

that allows pO)) to be ((32|) a Sturm Liouville well posed problem with weig ht function M^L = CQs\C^\pr/5) 
The ortonormalization condition to be satisfied by 0^ is 11] 



f 

Jo 



dr M^L 0;(r) 0,(r) = / dr cosh~^\/3r / 5) cf>*^ 0„ = 5^^, (35) 
10 Jo 

so the equivalent wavefunction in 1 dimensional quantum mechanics is 

ij^ir) = cosh-3*/2(/3r/,5) 0^(r). (36) 
Finally, the massless normalized equivalent wavefunction is given by: 

V>o(r) = koCoah-^^^^{l3r/6). (37) 
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This function is strongly decaying, as seen in FiglH So we conclude that the massless, spin two, gravitation mode is 
localized on the 3-brane and strongly concentrated around r = as we expected for a RS scenario. 

Although general wavefunction solutions are rather cumbersome, what is really important is the asymptotic behavior 
of the massive modes. Far from the vortex core or in the thin domain wall limit, when 5 — > 0, we can approximate 
(El by 

(38) 



So the massive wavefunction ([32)1 could be approximated by 

V;^(r)-e-i^'-</.^(r), 

and the localized zero mode by 



(39) 



(40) 



The massless zero mode is then localized in the vicinity of r = 0, that is on the 3 brane where the known universe 
is located, and decays exponentially when when r is increased. 

Using the approximation (|55)) and ([55)) into the differential equation and taking the limit /3r/S — > oo for which 
tanh(/?r/(5) -> 1, we obtain 



whose solution is given in term of bessel functions 



(41) 



(42) 



That eigenstates are not bounded to the brane and have infinite norm, in concordance with the result in 



FIG. 2: Graviton massless zero mode Ea. (|37p : /i = 0, = 1,5 = 1 



IV. MAXWELL EQUATIONS ON THE TOPOLOGICAL ABELIAN STRING 

We will now consider the background metric given by (El, where M{r) — L{r) is given by ([M|) . Performing variations 
in the gauge field, on the action ((T|) the field equations are: 



-G 



-GF 



AB 



+^l(0aS^_^gi3- 



(43) 



Here we will assume the Nielsen Olesen anzats ([6]) as valid only to zero order, and look for the equations of the first 
order fluctuations but with the simplification given by (I10[) . Moreover we will impose axial symmetry for all initial 
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fields, sources and potential, that means that both gravity and electromagnetic fields must have axial symmetry. Axial 
symmetry jointly with the abelian gauge transformation of the action implies we have two U{1), one for the gauge 
field and another for the invariance under spatial rotations in the 2 dimensional space with coordinates (r, 9) as in 
00. 



The gauge invariance allows us to impose two gauge fixings 







Ag ^ - 



(44) 
(45) 



as exact equations at all orders in perturbative theory. Equation is equivalent to pUj) . The 4 dimensional Maxwell 
field is considered to be a fiuctuation from its null background value ([6]) 

A^ = O + ^^(a;^r,0) 
III/ 4" d^Au di/A^ 

We will also assume (/) to be a background soliton, so we will look for solutions that preserve the background 
condition 



(j)^ v f{r) e*' 



(46) 



at least at first order. So the scalar field "current" contribution to the right term in (|43l) vanishes. 

Equation (|43p will be exact to first order in perturbative analysis for Ai_i with the metric given by the background 
([2]) and ((24|) . Note that in this approximation we have neglected graviphotons and graviscalars (spin 1 and from 
fluctuations of the metric) inspired by the results in [iB] and [l3| because graviphotons and graviscalars were found to 
be non normalizable either by their behavior at r — >■ or at r — >■ oo. The point of view here, is that the contribution of 
the non localized field to the electromagnetic fluctuation will be weak. As the norm of these wavefunctions divergesm, 
the value of the "normalized" wavefunction on the 3 brane and therefore its superposition with localized fields will 
be small. Although a renormalization procedure would exist, graviphoton and graviscalar are first order fiuctuation 
of the metric and as J^^^ is itself a first order, then the graviscalar and graviphotons interaction terms in ((43|) 
are second order terms and will not be considered here. 

The equation ((43l) for B = 9, due to (144]) and Q up to first order is: 



duJ'"'^ = 



a^a^^f = 



while the equation for B = r, due to (P5|) and ^ up to first order is: 



d'^dfiA'' = 



so both equations reduces to: 



d^A, 







(47) 



That is the 4 dimensional Lorenz gauge condition and comes from the fact that the electromagnetic U(l) gauge 
invariance was broken by the gauge fixings (j44l45p . 

If we add to the action ([1]) and interaction term between an external current coupled with the gauge field 



S^nt = J dx'^V^lABjgt] 



where this "external current" is given by 4 dimensional term "on shell" on the 3 brane: 
only the for the case B = fi equation (j43p will acquires a current term 



(48) 



(49) 



(50) 
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that could be break into two equations 



d^T^"" = (51) 

-1 



^^^dr[M'LdrA,] = j^d^eA.+v^e'WnfA^ (52) 

So we almost recover the 4 dimensional Maxwell equations with a 4 dimensional source on the brane, in the Lorenz 
gauge (|i7)) . if equation (15^ is accomplished. As the source in equation ([ST|) may depends on (r, 9) we are not exactly 
recovering Maxwell, unless the 4 dimensional current J^^ — J^jj{x^) depends only on or both J^^^^ and J^j^ has 
identical warping factors. Then we will recover completely the 4 dimensional Maxwell equations when r is close to 
zero. In this case the photon is simply the 4 dimensional vector potential multiplied by a warping factor (j60p that 
confines the photon to the 3 brane universe. To see that the former is indeed the case,we must proceed to the Fourier 
analysis of the equations in the following section. 

V. LOCALIZING THE PHOTON ZERO MODE ON THE TOPOLOGICAL ABELIAN STRING. 

If we expand the gauge field equations in a Fourier series as was done in a previous section and in ^] for the tensor 
case: 

^^(x^ r, 9) = A^{x) aiir) e'"" (53) 
then the equations for the Fourier coefficients using ([55)) in ([5^ is: 



1 d 



or 



qfai (54) 



with gf(r) = y^e'\\f{r)f-j^^^ (55) 

that is very similar to the equation (pO)) for the gravity case. 
Two important differences arise from equations (|54p and ([5(1)) : the first is that in ([50)) the square mass of the graviton 

= ml — (■^)^-|rj-isa constant eigenvalue due to M(r) — L{r) as in [llj, while for ([54| the "charge square" ([55)) is 

a complicate function. The second difference is that unlike the graviton case for which /z^ ^ 0, in equation ()54|) there 
is not physical reason to avoid g? < 0. In fact as lim L{r) = implies that asymptotically g? — oo for all non zero 

r— )-oc 

modes Z ^ 0. 

Equation ([54)) could be written for the massless case for all fouricr modes as 



M' 

a[' + 3— a[~qf ai^O (56) 

here we will take the values for M{r) from equation ()24l) and f{r) from ([22)) for the topological abelian string solution 
obtained in a previous section and in [23| : 

a'l' - 3l3ta.nh{pr/S)a[ - (^v^e^ ||/o arctan (sinh l3r/S)f '^^^^~Js)) " ^' ^^'^^ 

Equation ()57l) is rather involved even for the zero mode / = case, so we will use the same approximation ([38]) that 
was used in [11| to calculate the massive modes for the gravitational case, jointly with the approximation (see ([22])) 

/(r) ^ nVS (58) 

used in pol |. Both approximations are valid for r > where 5 is a not dimensional parameter related to the string 
thickness and /3 the Randall Sundrum warp factor. 
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First we will address the massless zero mode / = case, later in this section we will study the massless / ^ case. 
The massive case will be consider in the next section. Exact numerical solutions could be obtained as show Fig ([3]). 
For the / = massless case we have 



ap + 3/3 Op - (Jo ao = 0, 
that has the following base of solutions: 



ao(r) oc exp 



3/3 



± 



2 2 2 2 c- 

(7q = W e TT 0, 



(¥)• 



9o 



(59) 



(60) 



with the plus sign (j60p is diverging, with the minus sign is converging provided 



(5 < 



47r 



(61) 



it is also a non oscillating when (5 < j^-^r^^ that is fulfilled by the former equation ((6T|) that will be referred as the 
thin string limit. Note that in the thinnest limit, when (5 — > the zero mode converges to a constant value 

lim ao(r) — > 1 

but the coefficient can be normalized in the curved space we are working on. 
Equation (1541) could be rewritten as a Sturm Liouville equation: 



d_ 

dr 



+ q{r) ai = p{r) A a; 



(62) 



p{r)=M^L, q{r) = NrL v^e^\\f{r 



and 



p(r) = M(r) 



(63) 



where the eigenvalues A = P / are positive and real constants. The operator is self adjoint and the Sturm Lioville 
problem has solution with boundary conditions 



«;(o) = 



a;(oo)=0 



(64) 



that are of the same kind of condition taken in [Tl|. Note that the normalization factor comes from p{r) = M{r), the 
weight function, so the ortonormalization condition is: 



drM{r) ai(r) a„(r) = \\Ni\\^ Si^ 



\N,r - 



dr M{r) af{r) 



(65) 



Then the modes could be written in term of a equivalent normalized function in flat space 



VMM 



(66) 



That is consistent with the normalization used in 11] and [20] for the graviton case, but in graviton case the weight 
function was M'^L, instead of M{r). 

Maxwell equation (|5ip will be solved in term of the zero mode only, if we assume that (x" , r) has not dependence 
on the angular variable 9 



^^(x^r,0) = A^(x)a^(r) 



(67) 



where 
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a^(r) = ^ cosh-^/2(/3r/5), \\Nof = cosh-\pr/5) \\a,f , (68) 

here ao is given by (|60p in the thin hmit approximation (j6ip or could be solved numerically from (|57p. 
In the null thickness limit (5 — >■ we simply have 



lim a; 

■5^0 



N 



cxp[--rj 



(69) 



As the massless zero mode has a M^/^(r) warping factor, to obtain the equations ([ST|) . we must introduce the same 
warping factor to the current 



J,4(x^r) = J|'^(x^) Afi/2(r)5, 



(70) 



So, using (|67l) and ([70]) in ([51]) we recover completely the 4 dimensional Maxwell equations very close to the 3 brane. 
In this case the photon is simply the 4 dimensional vector potential multiplied by a warping factor (j60p that confines 
the photon to the 3 brane universe. 

To obtain the non zero Fourier modes, I 0, we could write equation ((56)) using (|38l) and (l58t as: 



3p o; + go a; = (2 e. ) ^ ai 



these equations has the following base of solutions in the thin limit approximation (j6ip : 

ai{r) cx i ei^'' / 



± 



2 2"! Br 



(71) 



(72) 



where the notation = In[x\ has been used for Bessel functions. 

These functions are not bounded for r — > cx), and could not be normalized because integrals in ([S5)) are diverging. 
They neither can be normalized using the graviton norm as in (|35p . So massless non zero modes are not 

localized on the world brane and wanders into the bulk. The massless zero mode is bounded to the 3-brane only if 
(j6ip is accomplished. Finally we want to stand out that if (j61l) is not accomplished, ao will be an exponential growing 
factor and the photon wavefunction will not longer be localized on the world brane, and scape into the bulk. 




FIG. 3: Numerical integration of 6a. (|57p for the massless modes: I = (left) and I = 1 (right) e = 0.1, d = 50, /? = 47r/3, fo = 1. 



VI. PROCA MASSIVE MODES AND THE CORRECTIONS TO THE COULOMB POTENTIAL 

The orthogonality of Fourier modes implies that equation (|50|) . and therefore (|5T|) and (|52p must be solved order 
by order. The complete solution will be a superposition of the eigenfunctions with eigenvalues ^ and Fourier index / 

4)(x^r,^)=A«(x)aKr)e^'^ (73) 
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that is slightly more general that ([53)) . and could be understood as a collection of vector potentials, one for each mode, 
with electromagnetic fields 



J■(')(x^r,0)=F«(x) aiir) e''^ 
Then equation ([50)1 with not current term ( Jg^-j — 0) could be written as 



1 d 



M2 



APL dr 



APL -[a,{r)] 



«V||/(r)|| 



Rl L^r) 



(74) 



a,(r)A«(a;) (75) 



the sign in — -^j^ is what leads to a non normalizable expressions as ([72 
Note that we can add to both right and left sides a Proca mass term: 



aiir) 



1 d 



A("(^) 



L2 



a,(r)A«(:r) (76) 



where ([58|) and Lorenz gauge (|47l) have been used. 

Former equation ((75)) could be split in similar way to (IFTI) and (15^ as 



□A(')+m2A(') = 0, 



1 9 
' M^L dr 



or 



ai{r). 



(77) 
(78) 



Equation ([77)1 is a Proca equation with 4 dimensional mass m for the photon. Equation ([75)1 gives a regularized 
version of (|54|) when — I^/Rq > and is a self adjoint Sturm Liouville equation. 
We could write equation ((5^ using [M] and approximations 



x;' - 3/3 a; + a; + (22*e2^'-)(m2 - j^) ai = 



these equations has the following base of solutions: 



„2 ^p/3r 



W(¥)- 



„2 ^p/3r 



(79) 



(80) 



where the eigenvalues /x —^m? — > are effective mass term and notation J[n,x\ = Jn[x\ and K[n, x] = Yn[x\ 

has been used for Bessel functions 

To obtain the first order electromagnetic field complete solution, we must sum over all the Fourier modes Z, and 
as ([77]) and ([75)1 depends on the mass, we must sum also sum over all the masses m > l/Ro, or equivalently on the 
eigenvalues fi as indicated: 



^.(x^r,0) = ^ J2 A«(:r)aKr)e''«^^^A«(x)a^(r) 
The expression for ([80]) could be simplified in null string thickness (5 — > to: 



iW 



(81) 



(82) 



where a more standard notation for Bessel functions has been introduced and and I?^ are constants. That 
expression is quite similar to that for gravitons [ll| and domain walls [1] [111 changing only the type of the Bessel 
function and the exponent in the warp factor. 

In the limit r — > oo, the solution (j82p grow exponentially for non zero fi. The standard way to regularize this, is to 
introduce a finite radial distance cutoff rmax where the boundary condition (|64l) will be imposed for r — rmax instead 
of r = oo 
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a;.(0)=0 a;,(w) = (83) 

These boundary conditions, jointly with the normalization condition (j65p leads to a discrete spectrum 



/x„ « /3(n + l)7r e-^'^— (84) 
for enough large n S Z, and solving for the constants we get: 

- ^-^i/^Iy] (85) 



When the cutoff is imposed the sum in Fourier expansion (1811) will be replaced by a sum over n, because (j84p are 
the only acceptable values of fi that will lead to renormalizable wave functions as expected by (j65l) . 
The normalization (I55t of the massive modes, can be calculated using ([5^ . ([55|) . as: 

||JV„f = I dr e^" (yy,[i^l J,/s[^^e* ] - J„,|^| Vi/^ly^* 

where the asymptotical form for ^3/2(^6'^'^] and Y3/2[^e^''], was used due to the large factor e^^ for large r — >■ rmax 
and enough large n 

When the cutoff is imposed, the sum in ((53)) and (|8T|) turns out to be over n e Z. So it will be used from now on 
in the rest of this work. Therefore, massive electromagnetic modes have fourier coefficients given by (|66p: 



<{r) 



£2^ 



(87) 



A Schrodinger equation could be obtained for the massive modes making the change: 

zi„(r) =a^(r) Af^ 
that accomplishes 

-ul,{r) + VEM-u,,{r) = (89) 

3M" 3 (M'V 2/ ^ 

The potential could be found calculating the integral of the Feynman graph at the tree level, for the interchange of 
a virtual photon between two stationary charges giand 52 placed on the 3-brane at r = 0, in the limit in which the 
photon energy A:° -> goes to zero: 



I n>0 J 



ViR) = ^ ll^^o(r = 0)||L.. +P/ > : > : hn(r = 0)|r^,,,exp(-m R) } (90) 



where R = y IjS'll is the spatial 3-distance between the charges, m = y/ /j,^ + P/Rg is the Proca mass, ||it|jg^^^ is the 
usual Schrodinger norm J dx'^u*u and pf is a polarization factor. 
For both massless and Proca photon the polarization factor is: 

lim p/ = g"" - = 1 (91) 



13 



where fc^ the momentum of the virtual photon, is the conjugate Fourier Transform variable to x^. 

Note that (|M| takes into account the contribution of the Maxwell or zero non massive propagator ^ and the 

contribution of the massive or Proca propagator "^^^"^ weighted by the squared norm of its wavefunction (j88p . 
that is its probability of been found at r = 0. Although the norm of the massive modes (|87l) is infinite, due to the 
large factor e^'""'"' in equation (f86|) . a finite limit could be found by considering the infinite sum in ([90|) .Using that 
lim[un(r)] — a^(0), because M(0) = 1, therefore the correction to the Coulomb potential could be calculated as: 

r— >0 

^ [ I n>o ||a^(r = 0)||g^^^ J 

This infinite sum could be transformed into an integral over the masses > 0, using the quantization condition (j84p 
for which: 

An = A/x„ ^ -Ig^'"— (93) 

So we get to an integral over the masses of the modes: 

|2 



The former integral could be obtained using (|86l) . ([57]) . ([55)1 and ([55]) . Taking either the large or small limit for 
in we get: 

llA^nir = ^ 2 2 (95) 

^ fJ-n 

Evaluating ([57)) at the brane, r = 0, and using (jHS)) and Bessel identities we simply obtain: 

l|a»(0)f. 



and using the thin limit for the zero mode ([69]) a^(0) = ^/^, we obtain: 



^^=^E/dM^ eM- ^/J^TpM R) (96) 
' 

Note that in the former integral the cutoff factor e^'"™"" was simplified, so the result is cutoff independent. Instead of 
equation (|69p , equations (ISUl) and (|68p could be used, making a Taylor expansion in S gives the same result at second 
order term in perturbative theory. In the sum there is implicit a step factor = 1, if we divide and multiply by Ro 
we get: 



AV 



gl 92 

R 



-y OO 

-1™ E§^ /dA* exp(-VM2 + P/i?2 i?) 



2Ro 
/?7r 



(97) 



That in the limit R^, >> A/ and making the change of variables 



Ro 



/.= X, 



Ro 
dfi = dY 



transform (|97p into a double integral 



Ay = 



gl g2 



^ J y dX dY exp(-VX2+Y2 i?) 





(98) 
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that could be calculated changing to polar coordinates with radius p = ^X? + and integrating over the area of 
the upper left quarter of the plane between [0,7r/2]. So finally: 



AV = 



R 



^ I dp cxp(-p R) 



gl g2 

R 



Ro 

13 i?2 



(99) 



If we have a non warped, flat space with 5 dimensional spatial coordinates, a would be expected for the Coulomb 
potential using Gauss law that is coincident with our result. 



VII. SUMMARY AND OUTLOOK 



In this work, the confinement of electromagnetic field is studied in axial symmetrical, 6D warped World Brane, using 
recently proposed [20| topological abelian string vortex solutions as background. The field equations were calculated 
only to first order in perturbative analysis for the vector field. The metric field was assumed to be exact and the 
solitonic scalar solution was assumed as solitonic background. As the theory is assumed to be axial symmetric and 
also containing classical electrodynamic or Maxwell theory, we have a U{1) x U{1) invariant theory as in [Tg'] and 
[l7j . that allow us to make two gauge U(l) fixings (|44l) and (|45|) . that is consistent with the topological abelian string 
solution, and lead to simplifications in the spin 1 fluctuation equations. 

There are several conclusions we found throughout this work: 1. There is a massless, spin 1, Fourier zero mode, 
bounded to the 3-brane universe, with a warping factor in the bulk if the string-vortex is thin enough. The shape of 
this mode could be seen in Fig([3]) from numerical integration, or obtained in the thin string approximation (|6ip using 
(p7)) . ([SO]) and ([55]) even could be simplified to ([S^ in the null thickness limit 5 — >■ 0. 

2. The massless zero mode is consistent with Maxwell equations (jSip with an external current warped throughout 
the bulk. 

3. All other modes in the photon expansion ([73| are not localized and massive. They follow Proca equations (|77p . 

4. The main conclusion of this work is that the correction to the Coulomb law, produced by the massive modes, in 
the thin string limit is: 



V{r) ^ ^ (l + ^) (100) 

The expected Coulomb potential for flat 6D spacetime is but in this case, it also depends on the factor that 
could be seen as the radius of compactification of the angular coordinate, and Rrs = ^ that is the distance at which 
the metric falls by a factor of exp(— 1) due to the warping factor in Randall Sundrum theories. The value of Rrs is 
unknown, and could be in a very wide range, as short as several Planck length up to the experimental limit for 
Newton potential around 10~^cm [l^. The value of Ro is also not know, but must be larger than Planck length, and 
shorter than the experimental limit of validity for the Coulomb law of 10~^^cm [2^ [2^. Experiments to test at short 
distances the electric Coulomb potential, could be used to establish the existence or not of extra warped dimensions 
at distances of 10~^^cm in the near future. That is an increase by a factor of 10^^ with respect to the actual capacity 
to observe gravitational effects due to the warping or compactness of extra dimensions. 

Observable effects of the correction to the Coulomb potential (|100p implies at least a deviation of = 0.1, 
because this ensures a new decimal figure that could not be explained by potential alone. In the most optimistic 
scenario: assuming Ro ~ lO^^^cm and ^y- = 0.1, in order to obtain an observable change in Coulomb law at a 
distance of i? « 10~^^cm (that could be achieve at the LHC) implies that Rrs must be greater than 10~^^cm. But if 
Ro ~ 10~^^cm and Rrs ~ lO^^'^cm, for a distance of i? « 10~^^cm then we obtain that = 10~^^ so the precision 
of the Coulomb law could be amazing. 

Proca photons has a long history in both theoretical and experimental physics, a striking limit of lO^'^^gr for the 
photon mass has been set by astronomical measures. As we have seen in the previous sections, most of the photons 
will correspond to the massless zero mode, that are localized at r = 0, just over the 3-brane universe. On the other 
hand, massive photons are not bounded to the 3-brane universe and wander in the extra dimensions, throughout the 
bulk, so the probability of catch one of the Proca photons is extremely low. Lab experiments [1^ and astronomical 
measurements [28*1 limiting the photon mass, assume that all photons has the same small mass. As this is not the 
case for this model, these mass limits do not apply. 
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